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The nutation (wobble) of a spinning spacecraft in the presence of energy dissipation is a well- 
known problem in dynamics and is of particular concern for space missions. The nutation of 
a spacecraft spinning about its minor axis typically grows exponentially and the rate of 
growth is characterized by the Nutation Time Constant (NTC). For launch vehicles using 
spin-stabilized upper stages, fuel slosh in the spacecraft propellant tanks is usually the 
primary source of energy dissipation. For analytical prediction of the NTC this fuel slosh is 
commonly modeled using simple mechanical analogies such as pendulums or rigid rotors 
coupled to the spacecraft. Identifying model parameter values which adequately represent 
the sloshing dynamics is the most important step in obtaining an accurate NTC estimate. 
Analytic determination of the slosh model parameters has met with mixed success and is 
made even more difficult by the introduction of propellant management devices and 
elastomeric diaphragms. By subjecting full-sized fuel tanks with actual flight fuel loads to 
motion similar to that experienced in flight and measuring the forces experienced by the 
tanks these parameters can be determined experimentally. Currently, the identification of 
the model parameters is a laborious trial-and -error process in which the equations of motion 
for the mechariical analog are hand-derived, evaluated, and their results are compared with 
the experimental results. The proposed research is an effort to automate the process of 
identifying the parameters of the slosh model using a MATLAB/SimMechanics-based 
computer simulation of the experimental setup. Different parameter estimation and 
optimization approaches are evaluated and compared in order to arrive at a reliable and 
effective parameter identification process. To evaluate each parameter identification 
approach, a simple one-degree-of-freedom pendulum experiment is constructed and motion 
is induced using an electric motor. By applying the estimation approach to a simple, 
accurately modeled system, its effectiveness and accuracy can be evaluated. The same 
experimental setup can then be used with fluid-filled tanks to further evaluate the 
effectiveness of the process. Ultimately, the proven process can be applied to the full-sized 
spinning experimental setup to quickly and accurately determine the slosh model 
parameters for a particular spacecraft mission. Automating the parameter identification 
process will save time, allow more changes to be made to proposed designs, and lower the 
cost in the initial design stages. 
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I. Introduction 


Spinning a spacecraft or an upper stage has been a long standing method for accomplishing stability and 
pointing accuracy of the spacecraft with a minimum of hardware, complexity and expense. While spinning a 
deployed spacecraft over its operational lifetime has generally fallen out of style in favor of the more modem three 
axis stabilized active systems popular today, there still is a community of users that have to deal with spin stabilized 
upper stage dynamics. Many NASA and DoD payloads are launched on Boeing Delta II expendable launch vehicles 
with a spinning solid rocket third stage. This particular version of the Delta II has been very popular for NASA 
interplanetary missions. Because of this, NASA’s Expendable Launch Vehicle program office at Kennedy Space 
Center has been investigating ways to improve their understanding and ability to model spinning upper stage 
dynamics. While this research work has important near term applications for expendable launch vehicles, it also has 
significant implications for NASA’s future manned space program. Spinning a large manned vehicle (or perhaps 
segments of one connected by a long tether) is the only practical way to obtain “artificial gravity”. Long duration 
space missions may require some form of artificial gravity to counteract the effects of extended weightlessness on 
the human body. 

Liquid slosh in the fuel tanks of an attached spacecraft has been a long standing concern for space missions 
with a spinning upper stage. Energy loss through the movement of the liquid fuel in the fuel tank affects the 
gyroscopic stability of the combined spacecraft and upper stage. Energy loss leads to an ever increasing wobble or 
“nutation” which can grow to cause severe control issues [1]. The more vigorous the slosh the greater the energy 
loss and hence the greater the nutation. The ’’nutation angle” is defined as the angular displacement between the 
principal axis of rotation of the spacecraft and its angular momentum vector and is a measurement of the magnitude 
of the nutation [12]. The amount of time it takes for the angle of the nutation to increase by a factor of e 1 is defined 
as the Nutation Time Constant (NTC) and is a primary factor in the determination of the long term stability of the 
spinning spacecraft during the upper stage burn. The NTC can sometimes be very tedious to calculate and very 
inaccurate if calculated during the early stages of spacecraft design. 

There is a high degree of uncertainty in predicting the effect of liquid fuel motion in spinning spacecraft. 
The resulting nutation growth can be excessive and can pose a threat to the mission. Missions have been lost 
because of excessive and unanticipated nutation growth (Explorer I, 1958 and ATS-5, 1969 being early examples). 
Purely analytical methods of predicting the influence of onboard liquids have been generally unsatisfactory [1]. The 
NTC values provided analytically are quite often significantly different than actual flight values. Hence, there is a 
need to identify conditions of resonance between nutational motion and liquid modes and to understand the general 
characteristics of the liquid motion that cause the problem in spinning spacecraft. The proposed research is a first 
step in trying to understand and model certain modes of induced resonance found during experimental testing and 
during flight. This study will focus on the modeling of fluid motion and will utilize the results obtained to develop a 
more accurate prediction of the fuel slosh effects on spin stabilization of spacecraft. 

During the initial design of spacecraft, use of purely analytical means of predicting the influence of 
onboard liquids has not worked well. Computational fluid dynamics software packages provide some insight, but it 
turns out that they have several shortcomings. Their complexity and inability to accurately model the coupling 
effects of sloshing mass on the six degree-of-freedom motion experienced by the spacecraft make their application 
problematic. 

Liquid oscillations in spinning tanks have been studied in the past. Liquid oscillations in spinning fuel tanks 
produce very different response characteristics compared to those of non-spinning fuel tanks [3]. An energy sink 
model was originally developed by Thomson [4] to include the effects of small, passive sources of energy 
dissipation. This model does not work well for spacecraft fuel slosh energy dissipation due to the fact that fuel mass 
is a large fraction of the total mass of the spacecraft. 

Extensive analysis has been done on the different tank shapes and locations, as well as the use of propellant 
management devices (PMD). A summary of this analysis, like that reported by Hubert [11] shows the vast 
differences in possible behaviors of different designs. For the off-spin-axis-mounted, cylindrical tanks with 
hemispherical end-caps that have been popular in a number of spacecraft programs, a number of relatively simple 
mechanical models have been developed. Hubert also notes that one of the most difficult aspects of employing such 
mechanical models is in the selection of appropriate parameters in the model. 

Use of mechanical analogs such as pendulums and rotors to simulate sloshing mass is a common alternative 
to fluid modeling. A homogeneous vortex model of liquid motions in spinning tanks and an equivalent mechanical 
rotor model was developed by Dodge et al., [5]. An approximate theory of oscillations that predicts the 
characteristics of the dominant inertial wave oscillation and the forces and moments on the tank are described. 
According to Dodge et al., the pendulum model simulates a motion that does not involve an oscillation of the center 
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of mass. Therefore, it is not a valid model of inertial wave oscillations. Weihs and Dodge [6] illustrate that the free 
surface effects can be ignored when the liquid depth is small. 

A 3-D pendulum model was proposed by Green et al. # [7]. There was evidence of liquid resonance from the 
experimental data. The resonance was closely tied to the tangential torque and to a lesser degree to the radial torque, 
and there was little or no resonance in the force measurements. Green at al., proposed a rotary oscillator concept to 
simulate the torque resonance in tangential and radial directions. This rotary oscillator model was superimposed. on 
the pendulum model to provide the overall response of liquid oscillation in the tank. 

The current research effort proposed is directed toward modeling the resonance in fuel slosh of spinning 
spacecraft using simple pendulum analogs. The pendulum analog will model a spherical tank with no PMD’s. An 
electric motor will induce the motion of the pendulum to simulate free surface slosh. Parameters describing the 
simple pendulum models will characterize the resonant modes of free surface sloshing motion. The one degree of 
freedom model will help to understand the fuel slosh resonance and serve as a stepping stone for future more 
complex simulations to predict the NTC accurately with less time and effort. 


II. Problem Definition 

Various simulation parameters are estimated by matching the pendulum/rotor model response to the 
experimental response of full sized test tanks in NASA’s Spinning Slosh Test Rig (SSTR) located at the Southwest 
Research Institute in San Antonio, Texas. The experimental set-up of the SSTR is shown in Figure 1. The SSTR can 
subject a test tank to a realistic nutation motion, in which the spin rate and the nutation frequency can be varied 
independently, with the spin rate chosen to create a centrifugal acceleration large enough to ensure that the 
configuration of the bladder (PMD) and liquid in the tank is nearly identical to the zero-g configuration. A complete 
description of the actual tests, data acquisition and analyses of data for the Contour mission is provided by Green, et 
al., [7]. The fuel motion is simulated using models with various parameters (inertia, springs, dampers, etc.,) and the 
problem reduces to a parameter estimation problem to match the experimental results obtained from the SSTR. 

The SSTR can accommodate a full-sized fuel tank complete with any internal PMD for testing. The SSTR 
measures and records the force and torque response of the fuel tank to the internal slosh motion of the fuel. It has 
the capability to identify and characterize slosh resonances. The data from the tests are used to derive model 
parameters that are then used in the slosh blocks of a MATLAB/SimMechanics-based spacecraft and upper stage 
simulation. Currently the identification of the model parameters is a laborious trial -and -error process in which the 
equations of motion for the mechanical analog are hand-derived, evaluated, and their results are compared with the 
experimental results. 

The proposed research is an effort to automate the process of slosh model parameter identification using a 
MATLAB/SimMechanics-based computer simulation of the experimental SSTR setup. Different parameter 
estimation and optimization approaches are being evaluated and compared in order to arrive at a reliable and 
effective parameter identification process. To evaluate each parameter identification approach, a simple one-degree- 
of-freedom pendulum experiment is being constructed and motion will be induced by an electric motor. By 
applying the estimation approach to a simple, accurately modeled system, its effectiveness and accuracy can be 
evaluated. The same experimental setup can then be used with fluid-filled tanks to further evaluate the effectiveness 
of the process. Ultimately, the proven process can be applied to the full sized SSTR setup to quickly and accurately 
determine the slosh model parameters for a particular spacecraft mission. 

The problem with modeling the complete SSTR as a starting point is that there is considerable complexity 
in the SSTR machine itself. By reducing the problem to that seen in Figure 2, a better understanding can be made of 
the effectiveness of the optimization and estimation approaches and to the fundamental slosh behaviors of the liquid 
without having to model all of the complexity of the SSTR first. The fixed weight mass represents the amount of 
fuel that is not undergoing free surface slosh while the free surface fuel slosh mass is represented by the mass 
attached to the pendulum. 
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Figure 1. Schematic diagram of Spinning Slosh Test Rig (SSTR) 
at Southwest Research Institute 



Figure 2. Schematic diagram of one degree of freedom (DOF) pendulum analog 

at Embry-Riddle Aeronautical University 
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III. Method of Approach 


A spherical tank with no PMD’s undergoing free surface slosh is the simplified model for the pendulum 
parameters. Free surface slosh occurs during small cyclic frequencies with a well defined resonant frequency [13]. 
The only sloshing motion assumed to be taking place in this simplified model is that of a surface wave. The rest of 
the liquid is at rest and can be treated like a fixed mass. Initial pendulum properties are found by the use of a 
program developed by Dodge [10]. This “Dodge” code predicts the modes of the fuel tank with that of a pendulum. 
The tank/fuel parameters such as shape and kinematic viscosity are provided as input to the program. Values for 
various pendulum parameters are provided on the output as seen in Figure 3. 


LIQUID MASS [mass units] 

LIQUID SURFACE HEIGHT above z=0 [length units] 
FIRST MODE PARAMETERS 

Pendulum mass [mass units] 

Pendulum length [length units] 

Pendulum hinge z-location [length units] 
Pendulum % critical damping 

Ratio of slosh amplitude to pendulum amplitude 
SECOND MODE PARAMETERS 

Pendulum mass [mass units] 

Pendulum length [length units] 

Pendulum hinge z-location [length units] 
Pendulum % critical damping 

Ratio of slosh amplitude to pendulum amplitude 
FIXED MASS PARAMETERS 
Mass [mass units] 

Z-location [length units] 

Mom. Inertia [mass*length A 2 units] 


= 7.400E+00 
= 1 . 524E-01 

= 4.289E+Q0 
= 9 . 771E-02 
= 1 . 527E-01 
= 5.030E-01 
= 1 . 340E+00 

= 1 . 075E-01 
= 2 . 888E-02 
= 1 . 483E-01 
= 5 . 030E-01 
= 3 . 281E-01 

= 3 . 004E+00 
= 1.519E-01 
= 2 . 791E-02 


Figure 3^ Example of “Dodge” code output for 12" sphere at a 50% fill level (water) 

units are in (m, kg, and s) 


The first mode parameters represent the majority of the fuel undergoing free surface slosh. The second 
mode parameters represent a small correction factor for the first mode parameters. This can be seen by observing 
that the fuel mass for the second mode parameter as an order of magnitude less than the first mode. This second 
mode mass can be added to the fixed mass with minimal concern for inducing error. After running the code for 
several fill levels, several plots can be created from the data. Figure 4 shows the weight distribution for various fill 
level percentages for a 12 inch sphere with water as the liquid. The code predicts that the maximum sloshing mass 
will occur at approximately 60% fill level. Figure 5 indicates the various pendulum lengths that are required for 
different fill levels. 
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Figure 4. Liquid mass distribution for 12" sphere 
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Figure 5. Pendulum length distribution for 12” sphere 
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Using the code’s data distributions along with the geometric/material characteristics obtained from the 
experimental setup (Figure 2), a computer simulation of the one DOF pendulum analog can be developed using 
SimMechanics software [8, 9] as illustrated in Figure 6. Each of the different parts of the simulation model in 
Figure 8 is located in one of the following four groups: 

Group I: 

This group simulated the electric motor and locomotive arm assembly. The locative arm consists 
of five different parts starting with the DC motor. These are the flywheel, flywheel linkage, 
locomotive arm (piston), and the “stinger”. The “stinger” is a flexible metal rod designed to 
absorb forces that are not coincident with the axis of the locomotive arm. Geometric parameters 
such as component mass and moments of inertia are fixed in this group. Desired Frequencies and 
simulation running time are also input parameters for this group. 


"Stinger" 


To Pendulum Frame 



Figure 6. Group I assembly 



Figure 7. Groups I-IV locations 


Group II: 

The location of the force transducer in Figure 7 is represented by this group. The force transducer 
connects the pendulum frame to the locomotive arm assembly. The total forces of the fixed and 
sloshing masses are recorded in this group. 
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Group III: 

Parameters in this group include non-sloshing fuel mass assumed to be fixed along with its 
geometric properties of the frame assembly illustrated in Figure 7. The frame is constructed out of 
1/8” aluminum and weighs of 1.46 kg. 

Group IV: 

Group IV simulates the sloshing fuel and is considered to be critical in the parameter estimation 
process. Fuel mass, hinge damping, and pendulum length are all parameters in this group. 
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Figure 8. SimMechanics model of a pendulum as seen from Figures 2 & 7 

As a proof-of-concept, results from this SimMechanics simulation incorporating the one DOF pendulum 
analog matches those predicted by the Dodge code for a 60% fill level as shown in Figures 9 and 10. Figure 9 is a 
force simulation that would be experienced by the force transducer in the experiment. The simulated locomotive 
assembly (Group I) is driven by parameters determined by the user. At the start of the simulation (t=0s) , the 
frequency is 0.5 hertz (30 RPM) and at time equal to 500 seconds, the frequency is equal to 3.0 hertz (180 RPM). 
At approximately 250 seconds, the pendulum reaches its first mode at a frequency of approximately 1.71 hertz. This 
frequency matches the natural frequency prediction of the “Dodge” code. The natural frequency is inversely 
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proportional to the square root of the length of the pendulum. . For the case of a 60% fill level, the length is 
85.05mm. 



Figure 9. SimMechanics output for sphere at 60% fill level 
force (lb) vs. time (s) 



Restore saved axes setting 


Figure 10. SimMechanics output for sphere at 60% fill level 
Rotational speed (RPM) vs. time (s) 

After confidence is gained with a 1 degree of freedom pendulum model, a spherical tank with fuel 
properties which match the Dodge parameters (water in this case) will then be tested experimentally to verify the 
validity of the pendulum analog as shown in Figure 11. The block diagram of the parameter estimation process is 
illustrated in Figure 12. 
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Figure 11. Schematic diagram of one DOF fuel analog 
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Figure 12. Block diagram of the parameter estimation process 
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IV. Conclusion 


The effects of fuel slosh aboard spinning spacecraft need to be accurately predicted to avoid mission 
failures. Using a combination of test derived fuel slosh parameters and computer simulations of the spacecraft 
dynamics, an improvement in the current ability to make predictions of NTC can be achieved. By applying the 
parameter estimation approach to a simple, accurately modeled system during initial stages of design, its 
effectiveness and accuracy can be evaluated. The same experimental setup can then be used with fluid-filled tanks 
to further evaluate the effectiveness of the process. Ultimately, the proven process can be applied to the full sized 
spinning experimental setup to quickly and accurately determine the slosh model parameters for a particular 
spacecraft mission. Automating the parameter identification process will save time, allow more changes to be made 
to proposed designs, and lower the cost in the initial design stages. Applications of an automated process to find the 
NTC will benefit all space exploration missions involving a spinning spacecraft. On future manned space 
exploration missions involving artificial gravity, understanding and being able to confidently predict the stability of 
the spinning human habitat will be crucial for the success of the mission. 
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